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Abstract 

In this paper we study the characteristic or generating function of a certain dis- 
continuous linear statistics of the Laguerre unitary ensembles and show that this is a 
particular fifth Painleve transcendant in the variable t, the position of the discontinu- 
ity. 

The proof of the ladder operators adapted to orthogonal polynomial with discontin- 
uous weight announced sometime ago [T3] is presented here, followed by the non-linear 
difference equations satisfied by two auxiliary quantities and the derivation of the 
Painleve equation. 
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1 Introduction 



In the theory of random matrix ensembles with unitary symmetry the real eigenvalues 
{xj}^=i have the joint probability distribution 



1 " 

P{Xi,...,Xn)dXi...dXn = , J]^ {Xk - Xj)^ Yi^oi^l)^^^ 



n\Dn 



:i.ii 



l<j<k<n 



1=1 



where Wq{x) with x G [a, 6] say, is strictly positive and satisfies a Lipshitz condition and 
has finite moments, that is, the existence of the integrals. 



x^WQ{x)dx, j G {0,1,2,...}. 
Here D„ is the normalization constant 



so that 



If " 

Dn[wo] = —; TT {xk-XjYY\^o{xi)dxi, 

^' >^K^]" l<,<fc<n ;=1 

/ P{xi, ...,Xn)dXi...dXn = I- 

J[a,b]" 



(1.2) 



:i.3) 



We include the cases where a may be — oo and/or b may be cxd . For a comprehensive 
study of the theory of random matrices see [2D] . 

A linear statistics is a linear sum of a certain function g of the random variable Xj : 



^9{xk). 



:i.4) 



fc=i 



The generating function of such a linear statistics is, by definition the average of exp(A X^fc fi'(^fe))' 
with respect to the joint probability distribution (1.1), where A is a parameter, reads. 



exp 



a,b]" 



More generally, we can consider 



^Yl Si^k) 



k=l 



P{Xi,...,Xn)WdXk. 



:i.5) 



fc=l 



/ 



WfM 



k=l 



P{Xi,...,Xn)WdXk 



:i.6) 



fc=l 



If / is a "smooth" function then asymptotic formulas for large n for the characteristic 
functions have been obtained for the Hermite case, where wq{x) = e~^ , x G (—00,00) 



by Kac [16] and AkhiezerfTj and generalized by many authors. See [8] for a history of 
this problem. These results are continuous analogs of the classical Szego limit theorem on 
Toeplitz determinants 

In the Laguerre case where wq{x) = x"e~^, x G [0, oo) an analogous formula was found 
recently for smooth functions in [4J. However, results for the situations where / has dis- 
continuities are harder to come by. We mention here the original studies in [6] where / has 
several discontinuities and which corresponds to the Hermite case. More general results can 
be found in P [101 El and [2T]. Also, results that correspond to a = ±1/2 and large n 
appear in [5]. 

In this paper we investigate the case where n is finite, and / is constant except for a 
jump at t G [0, oo), and is of the form 



/(x,t) =A + Be{x-t) 



:i.7) 



where 9{x) is one for a; > and zero otherwise and A > and -B > . In the special case 
of linear statistics the function g will take the form 



g{x,t):=e{x-t) ln(l + ^ 



0(t - x) In ( 1 - - 



1.8) 



where — 1 < | < 1. This corresponds to a function / where 



/(x,t) 



1- 



/3 



that is 



and 



+ 



A 



1 + 



P 



e{x-t) 



B 



-fr-(-r^ 



We also point out that if A = and B = 1 then we have the important case where we 
are computing the probability that all the eigenvalues are in the interval [t, oo). This case 
of course is what leads to the now well-known Tracy- Widom laws. More will be said about 
this later. 

Our main tool will be to use the theory of orthogonal polynomials. Previously, in random 
matrix theory one made use of the orthogonal polynomials associated to the weight that 
defined the ensemble. Fundamental quantities were then described in terms of Fredholm 
determinants. While both the authors are very fond of determinants, in this work, we do not 
consider Fredholm determinants. Instead we consider the polynomials that are orthogonal to 
the perturbed weight, that is a regular or "nice" weight multiplied by the discontinuous factor 
given in (1.7). In this manner we are able to use the results of the orthogonal polynomials 
to derive equations associated with the various statistics of interest. 



The idea is that we write the multiple integral in (1.6) as a Hankel determinant. We then 
need to know information about the norms of the orthogonal polynomials. To understand 
this we need to know something about the recursion coefficients of the polynomials. This 
will lead us naturally to another pair of auxiliary quantities that depend on t and n . In the 
paper they are called r„ and Rn ■ Using these auxiliary quantities we are able to produce 
the second order non-linear differential equations satisfied by S'„ = 1 — 1/-R„ which turns 
out to be a particular fifth Painleve transcendent, in addition to the Jimo-Miwa-Okamoto 
a form [TB] satisfied by the logarithmic derivative of the Hankel determinant with respect 
to t. We also wish to emphasize that the logarithmic derivative of the Hankel determinant 
can be computed very naturally in terms of our quantity r„(t) and its derivative and the 
relationships between these quantities arise naturally using this approach. 

We also derive a discrete version of the a form of a nonlinear second order difference 
equation satisfied by the same logarithmic derivative. Our computations show that in fact 
the values of our generalized polynomials at the end points of the intervals are intimately 
related to the resolvent kernels found in the standard approach of Tracy and Widom. This 
is really not surprising, since we are all starting with the same multiple integral. Rather, our 
point is that computations can all be made by using only the very basic theory of orthogonal 
polynomials. 

The Palinleve equation can be found in [2Z]- The second order difference equation [(4.29), 
Theorem 8], as far as we know is a new equation. 

In the next section the proof for a pair of ladder operators, and the associated supple- 
mentary conditions adapted to orthogonal polynomials with discontinuous weights which 
was announced sometime ago [13] will be provided. In section 3, a system of difference equa- 
tions satisfied by two auxiliary quantities r„ and i?„ (these will ultimately determine the 
recurrence coefficients for the orthogonal polynomials) are derived. In section 4 we derive a 
second order non-linear differential equation which turns out to be a particular fifth Painleve 
transcendent. In the process we identity that the quantity 

Sn{t) := 1 



Rnity 

to be such an equation. Furthermore we show that the logarithmic derivative of the gener- 
ating function 

Hn{t) := t- In G(n, t) = t- In D„(t) 
at at 

satisfies both the continuous and discrete a form of Painleve V. 



2 Ladder operators and supplementary conditions 

According the general theory of orthogonal polynomials of one variable, for a generic weight 
w , the normalization constant (1.2) has the two more alternative representations 

Dn[w] := det(/ii+j)^-io := det I / x'^^ w{x)dx\ (2.1) 

\Ja J jj=0 

n-1 

= n ^^■' (2.2) 

where the determinant of the moment matrix (/ij+j) is the Hankel determinant. Here 
{/ij}"=o is ^he square of the L^ norm of the sequence of (monic-)polynomials {-Pj(a:)}"=o 
orthogonal with respect to w over [a, 6]; 

h 
Pi{x)Pj{x)w{x)dx = Sijhj. (2.3) 

Therefore with reference to (1.2) and (1.5) the quantity that we need to compute is 

Dn[w] _mjh,it) 



G(t,n) 



^nH UZ'h. ' 



li=0 '"« 

where w{x,t) := x°'e~^{A + B9{x — t)) and hk(t) is defined by 

{Pkix)y{A + Be{x- t)) x" e-^'dx = hk{t). (2.4) 

We also denote 

Dr,{t) := Dn[w{.,t)\. 

This leads to the generic problem of the characterization of polynomials orthogonal with 
respect to "smooth" weights ^0(2^) perturbed by a jump factor where the discontinuity is 
at t. So if we write 

wj{x,t):= A + Be{x-t), A>0, A + B>Q (2.5) 

then 

b 

Pi{x)Pj{x)wo{x)wj{x,t)dx = 6ijhj(t). (2.6) 

It follows from the orthogonality relations that, 

zPn{z) = P^+,{z) + an{t)P^{z) + /?„(t)P„_i (2;). (2.7) 



This three term recurrence relations, together with the "initial" conditions, 

Po(^) = l, /5oP-i(;2)=0, 
generates the monic polynomials, 

P„(^) = ^" + Pi(n,t)^"-^ + ... (2.8) 

the first two of which are 

Po{z) = 1 

P,{z) = z-ao{t) = z-^^. (2.9) 

Note that due to the t dependence of the weight, the coefficients of the polynomials and 
the recurrence coefficients a^ and /?„ also depend on t the position of the jump. However, 
unless it is required we do not display the t dependence. 
From (2.7) and (2.8), we find, for n e {0, 1, 2, ..} 

a„ = pi(ra,t) - pi(ra + l,t), 

n-l 

J]«, = -Pi(n,t) (2.10) 

i=o 

where Pi(0,t) := 0. 

From (2.6) and (2.7) we have the well-known strictly positive expression, 

/?n:=T^. (2.11) 

n-n-i 

Another consequence of the recurrence relation is the Christoffel-Darboux formula 

^ hk K-i{x-y) 

The above basic information about orthogonal polynomials can be found in [26] . 

In this section, we give an account of a recursive algorithm for the determination of the 
an, Pn for a given weight. This is based on a pair of ladder operators and the associated 
supplementary conditions to be denoted as (^i) and (^2) . For an general "smooth" weight 
the lowering and raising operators has been derived by many authors [Sj [TTl [121 l2l]- We 
should like to note here A. P. Magnus's contribution to this formalism [TTl [13 [19] . Indeed, 
we have been motivated by the investigation of [i9l where he obtained the large n behavior 
of the recurrence coefficients of a generalization of the Jacobi polynomials in which the 
standard Jacobi weight is perturbed by a "line" analogue to the Fisher-Hartwig singularity. 
We end the discussion about the ladder operators with the remark that the supplementary 
conditions for orthogonal polynomials on the unit circle was found in [2] and have been used 
to compute explicitly the Toeplitz determinants with Fisher-Hartwig symbols. 

The lemma below gives a detailed proof of the ladder operators in the discontinuous case 
where the results were announced sometime ago |13j . 



Lemma 1 Let Wo{x), x G [a, b] be a smooth weight function where the associated moments, 

x^Wo{x)dx, J €{0,1,2,..} (2.12) 



of all order exist. 

Let Wo{a) = Wo{b) = 0, and vq{x) := —\nwo{x). 

The lowering and raising operators for polynomials orthogonal with respect to 

w{x) ■.= Wo{x)wj{x,t), 



are 



where 



P'^z) = -B^{z)Pn{z)+P^A^{z)Pn-i{z), (2.13) 

P'^_^{z) = [Bn{z) + v'^{z)]Pn-i{z)-An^i{z)Pn{z), (2.14) 



An{z) := -^ + 7- -^^ '-^P^{yw{y)dy 2.15 

z-t hn Ja z-y 

Bn{z) := ^ + ^ r "^'^"^ ~ '''^y^ Pn{y)Pn-Mw{y)dy (2.16) 

Z - t hn-l J a Z-y 

Rnit) := B^{Pn{t,t)r (2.17) 

r„(t) := 5-^^P„(t,t)P„_i(t,t). (2.18) 



where 



P„(t,t):=P„(2;,t)|^^^. 
Here \nwo{x), is well defined since Wo{x) is suppose to be strictly positive for x G [a, b]. 

Proof: We start from 

n— 1 
PL{z) = J2Cr.kPk{z), 
fc=0 

where Cnk is determined from the orthogonahty relations, 

r-6 



Cnk = ^ [ P'MPk{y)w{y)dy. 

hk J a 



Therefore 



ra-l 



PUz) 



^h Pniy)Pk{y)w{y)dy 



n-l ^h 

E 

k=0 
b 



Pk{z) 

a hk 



PMiPkivMy) + Pkiy) [B6iy - t)w^{y) + w'^{y)wj{y, t)\}dy 



n-l 

Pn{y) J2 

k=0 

n-l 



Pk{z)Pk{y) 

hk 



Bwo{y)S{y-t) + ^w{y) 



dy 



Pn{y) Yl ^^^^^{B wo{y)5iy - t) + K(^) - v',{y)]w{y)}dy 



fc=0 



hi 



Pn{y) 



Pn{z)Pn-l{y) - Pn{y)Pn-l{z) 



{B5{y - t)wo{y) + Kiz) - yM]w{y)}dy 



hn-i{z -y) 
where we have used integration by parts, (C-D), the definition of vq, (2.11) and that 

Pniy)Pkiy)wiy)dy = 0, A; = 0, 1, 2, ..., n - 1, 



to arrive at the above. A httle simphfication produces (2.15) and (2.16) follows from straight 
forward application of the recurrence relations. □ 

Remark 1. If Wo{a) ^ 0, WQ{b) ^ 0, the terms 



w{y) 



{Pn{y,t)y 



hn{t){z-y) 



and w{y) 



Pniy,t)Pn-iiy,t) 



y=a 



hn-i{t){z-y) 



y=a 



are to be added into the definition of An{z) and Bn{z) respectively. 

Remark 2. If there are several jumps at ti, ...,tN then the first term of (2.15) and (2.16) 

should be replaced by 



N 



E 



^n,j\S'j^ 'jj 



N 

^ z-t. 

j=l 3 



'njy^ji ''J 



where 



'°°'^*'kp,(ty,t)V 






Kit) 

Wo{tj 



'K-i{t)' 

t := (ti, ..., tAr). 



As in the case of the smooth weight the "coefficients" An{z) and Bn{z) that appear in 
the ladder operators satisfy two identities vahd for all z E CU {oo}, which we gather in the 
next lemma. 

Lemma 2 The functions An{z) and Bn{z) satisfy the following identities which hold for 
all z : 

Bn+i{z) + Bn{z) = {z- an)An{z) - \/q{z) {Si) 

l + {z- an)[Bn+i{z) - Bn{z)] = /?„+i A„+i (2;) - /3„A„_i(2;) (5*2) 

Proof: By a direct computation using the definition of An{z) and Bn{z) . □ 

It turns out that a suitable combination of (5*1) and (5*2) produces an identity involving 
^""q Aj{z), from which further insight into the recurrence coefficients may be gained. 

Lemma 3 An{z), Bn{z) and Yll=o ^ji^) satisfy the identity 

n-l 
[Bn{z)]^ + ^/q{z) Bn{z) +J2^ji^) = f3nAn{z)An-liz) (S!,) 

j=0 

Proof: Multiply (5*2) by An{z) and replace {z — q;„)A„(z) in the resulting equation by 
Bn+i{z) + Bn{z) + y'oiz). See (5*1) . This leads to 

[Bn+l{z)]^ - [Bn{z)]^ + Vo{z)[Bn+l{z) - Bn{z)] + An{z) = f3n+lAn+l{z)An{z) - (3nAn{z) An-l{z) 

Taking a telescopic sum of the above equation from to n — 1 with the "initial" conditions, 
Bq{z) = and /3oA_i{z) = 0, we have {S'2). □ 

Let y = Pn{z) we find by eliminating P„_i(z) from the raising and lowering operators, 
the second order differential equation 



Lemma 4 



y'\z) - ( ^(.) + ^ 



^) y'i^) + (b'^{^) - ^"^i^ + E^. W j yi^) = 0- (2-19) 



Proof: By a straight forward computation using (2.15), (2.16) and (S'2) . □ 

Recalling (2.17) and (2.18) we note that if Vq(2;) is rational in z then the difference 

kernel, [vq^z) — y'Q{y)]/{z — y) is rational in z and y. Consequently (5*1) and (S'2) maybe 

put to good use to obtain a system of difference equations satisfied by the auxiliary quantities 

Rn and r„ and the recurrence coefficients q;„ and /3„. This will be clear in the next section. 



3 Recurrence coefficients and difference equations. 

For the problem at hand, 

Wo{x) = x°'e~^, X G [0, oo), 

Vo(a;) := — \nwo{x) = —alnx + x 
and for a > 0, wo{0) = 0. Note that w{oo) = 0. An easy computation gives, 

Vo(^)-Vo(l/) ^ o^ 
z-y zy' 

Using these and integration by parts we have the following 

Lemma 5 

A^{z) = Ml + l^3M^ (3.1) 

z — t z 

where 
i?„(t) := Bw,{t)^^^^ (3.3) 

r„(t) := Bw,{t f-^^^l^^--;}^^^\ (3.4) 

hn-l[t) 

Proof: Through integration by parts we find. 



oo 



a i/"-V^ti;j(i/;t){P„(t/,t)}% = hn{t) - B w,{t){P^{t,t)Y (3.5) 

/•oo 

a I y"'^e'ywj{y;t)Pn{y,t)Pn-i{y,t)dy = -n/i„_i(t) 
Jo 

- Bwo{t)Pn{t,t)Pn-l{t,t), (3.6) 

and we have used the fact that 

d 
—Pn{x, t) = n Pn-i{x, t) + lower degree 

to arrived at (3.6). From (3.5) and (3.6) and the definitions of An{z) and Bn{z), (3.1) — (3.4) 
follows. □ 

Substituting (3.1) and (3.2) into (5*1) we find by equating the residues 

r„+i + r„ = Rn(t- an) (3.7) 

-(r„+i + r„) = 2n+ l + a-a„(l -i?„). (3.8) 



10 



Lemma 6 



^n+1 + ^ra 



2n+l + a + tR„ 

Rn{t - an) 



(3.9) 
(3.10) 



Proof: (3.7) + (3.8) implies (3.9) and we restate (3.7) as (3.10). □ 

Substituting (3.1) and (3.2) into (5*2), we find, after some elementary but messy compu- 
tations, 



n— 1 



[Bn{z)f + y'^{z)Bn{z) + 5^ A, 



rl {n + rn){a + n + rn) 



j=0 



(z - ty 






1 _ a _ 2(n+rn) 
t t 



z-t 



+ 



n 



-5^i?, + (n + r„)(^^-l^ + 



avr. 



and 



(3nAn{z)An^l{z) 



+ 



PnRnRn-l /3n(l — -Rn)(l ~ Rn-1 

{z-tf ^2 

1/1 1 



t \z — t z 
Now (5*2) implies 
Lemma 7 For a fixed t, the quantities r„, i?„, /5„ satisfy the equations 

^n Pn-f^n-'^ra— 1 



(n + r„)(n + a + r„) = /5„(1 - i?„)(l - i?„_i) 



n— 1 
j=0 



■J T" ^n 



a _ 2{n + rr. 
t ~ t 



(3.11) 



Pn [(1 - i?„)i?n-l + (1 - Rn-l)Rn] ■ (3.12) 



(3.13) 
(3.14) 



^ [(1 - i?„)i?„-i + (1 - i?„_i)i?„] . (3.15) 



Proof: The equations (3.13) — (3.15) are obtained by equating residues of (5*2). □ 

In the next Lemma an expression is found for /3„ in terms of r„ and /?„. 



Lemma 8 In terms of Vn and Rn, /3„ the off-diagonal recurrence coefficient reads 



Pn 



1 — R„ 



r„(2n + a)+n{n + a) + -^ 

Rn 



(3.16) 



11 



Proof: We eliminate /SnRnRn-i from (3.13) and (3.14) to find 

rn{2n + a) + n{n + a) = /9„(1 — /?„ — 

= /9„(1 



Rfi-i, 



(3.17) 
(3.18) 



D 



Rn) - -^■ 

In the last step we have used (3.13) to replace PnRn-i by r^/Rn- 

We note that -B > can always be satisfied for the proper range of A. 
The equation (3.9) states that «„ is linear in i?„ up to a linear form in n , together 
with (3.10) and (3.16), when combined with say, (3.13) provide us with a pair of non-linear 
difference equations satisfied by r„ and i?„. We state this in the next theorem. 

Theorem 1 The quantities r„ and Rn satisfy the difference equations; 

Rn{t — 2n — a — 1 — tR„ 



with the "initial" conditions, 



1 1 

R, 



R 



-n— 1 



r„(2n + a) + n{n + a) 



Ro{t) 
hoit) 





5 re-* 



1 



r(l + a) 



1 



X e 



(3.19) 
(3.20) 

(3.21) 
(3.22) 

'dx. (3.23) 

D 



Proof: This is simply a restatement of (3.10) and (3.13) with (3.9) and (3.16) 

We shall see that {S'2) automatically performs finite sums in "local" form, of the quan- 
tities Rn and «„. This will be seen later to be relevant in the evaluation of the derivative 
of \\iDn{t) with respect to t and the derivation of the Painleve transcendent. 

Theorem 2 



n-l 



i=o 



n-l 
i=0 



-t Tn — n{n + a) + P„ 

(3n 



-PiW 



Z Tn' 



(3.24) 
(3.25) 



Proof: From (3.15) we have. 



n— 1 

tJ^Rj 

j=0 



Tn [2(n + rn)+a-t]+Pn [Rn + Rn-1 ' 2i?„i?, 



n-'^n— 11 



r„[2(n + r„) 
r„[2(n + r„) 



a 
a 



t] + (3n [Rn + Rn-l] - 2< 

t] + j3n - rn{2n + a) - n{n + a) - 2rl 



-tTr. 



n{n + a) + I3n 



(3.26) 



12 



The second equality of (3.26) follows from (3.13) and the third equality follows from (3.17). 
The equation (3.25) follows from (3.9) and the second equality of (2.10). □ 



4 Py(0,-^,2n+l+a,-i) 



In this section we shall discover which of the auxiliary quantities defined as the residues of 
the rational functions An{z) and Bn{z) is a Painleve transcendent. 

This will be obtained from a pair of Toda equations which shows that the Hankel deter- 
minant is the r— function and these when suitably combined with the difference equations 
produce our Py- 

Taking the derivative of /i„(t) with respect to t, we find 

|ln/^„(t) = -Bw,{t)^^^^ = -R^{t), (4.1) 

and consequently we have the Theorem 
Theorem 3 

-t-\nD^{t) = -tJ2 In h,{t) 

j=0 
n-l 

= t'^Rj = -Pi{n,t) -n{n + a). (4.2) 

j=o 

Proof: The proof is obvious. □ 

The next lemma gives the derivative of Pi{n,t) with respect to t. 



Lemma 9 



-Pi(n,t)=r„(t). (4.3) 



Proof: Note the t dependence of p^ (n, t) . Taking a derivative of 

POO 

0= / Pn{x)Pn~l{x)wj{x,t)Wo{x)dx, 

Jo 
with respect to t, produces, 



= -B W^{t) Pn{t,t)Pr,.^{t,t) + f 

Jo 



-Pi(n,t)x"-i + 



Pn-l{x)Wj{x, t)wo{x)dx 



-B Wo{t)Pn{t, t)Pn-l{t, t) + hn-l J^Pli^, t) . 



and (4.3) follows. □ 

We expect Dnit) to satisfied the Toda molecule equation [25] and this should indicate 
the emergence of a Painleve transcendant. The question that we will address is "Which 
quantity is satisfied by this particular Painleve transcendant?" 

13 



Theorem 4 The Hankel determinant Dn{t) satisfy the following differential- difference or 
the Toda molecule equation f2^, 



,2d\ ^ ,.. / ^ ^ ^ Dn+l{t)Dn-l{t) 

t -In D4t) = -n{n + a) + ^^^ . 

Proof: Taking a derivative of (4.2) with respect to t and (4.3) imply 

Now substitute r„ given above into (3.24) to find, 



(4.4) 



r„. 



n— 1 , 

a 



t±lnDJt) 



^^ ^ dt 

j=Q 

= -t-\xiDn{t) 



- n{n + a) + Pn- 



dt 
where the last equality comes from (4.2). The equation (4.4) follows if we recall 

Pn 



hn-1 D 



since Dn = ho...hn-i- □ 

We now state a pair of somewhat non-standard Toda equations. 

Lemma 10 The recurrence coefficients an and (3n satisfy for n G {1,2,..} the differential- 
difference equations 

/3',it) = {Rn-l - Rn)Pn (Ti) 

a'nit) = rn- rn-1, (T2) 

with ro(t) and Ro{t) given by (3.22) and (3.23) respectively. 

Proof: These equations are an immediate consequence of (4.1), (2.11), (4.3) and the first 
equality (2.10). □ 

To discover the Py of our problem. We first state two preliminary lemmas describing 
the t evolution of r„ and i?„. 

Lemma 11 For a fixed n, Rnif) satisfies the Riccati equation, 

tR'^ = 2rn+{2n + a-t + tRn)Rn. (4.5) 

Proof: We begin with (T2) and replace r„+i by Rnit — an) — r„. See (3.7). This leaves 

a'^ = 2r„ - (t - an)Rn 

After eliminating a„ in favor of Rn with (3.9) we have (4.5). □ 
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Lemma 12 For a fixed n, r„(t) satisfy the Riccati equation, 

K = r}~ ^^: , {r^r - (2n + a)^ - n{n + a)-\. (4.6) 

Proof: By equating (3.24) to the last equality of (4.2), we find 

Pi(n,t) =trn- (3n- 
Taking a derivative of the above equation with respect to t and noting (4.3) we see that 

K = P'n 

= [Rn-1 — RnlPn 

tin 

and use have been made of (T2) and (3.13) to obtain the last two equalities. The equation 

(4.6) follows if we express /5„ in terms of r„ and Rn using (3.16). □ 

The next theorem shows that R^ is up to a linear fractional transformation a particular 

Py. 



Theorem 5 The quantity 



S„W:=l-^, (4.7) 

satisfies 

S'^ = ^'^""^ (S'J"^ -^- ^ {Sn-lf ^^2n + l + a)—-- ^^^^r, + I) ^ ^^ g^^ 

2>-'n(,-'- ^n) t Z t On t Z O^ J- 

wMch IS Py(0, -aV2, 2n + 1 + a, -1/2). 

In terms of the recurrence coefficient anif), we have, 

^' an{t) - {2n + a + 1) ^ ' 

Proof: Eliminate r„(t) from (4.5) and (4.6) and with i?„ = 1/(1 — Sn) gives (4.8). We 
have followed the convention of [14j. □ 

Remark 3. Note that for n = 0, (4.8) is satisfied by 

1 



So{t) 



Roity 

where -Ro(^) is given by (3.22) and (3.23) and ultimately in terms of an Incomplete Gamma 
function — a special case of the Kummer function of the second kind. Furthermore, since 
ro(t) = 0, it can be verified that -Ro(^) also satisfy (4.5) at n = 0. 

We may express the logarithmic derivative of Dn{t) with respect to t in the so-called 
Jimbo-Miwa-Okamoto a form. This is described in the next theorem. 
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Theorem 6 Let 

H„{t):=tj^\nD^{t), (4.10) 

then 

{tH'^f = A{Hlf[Hn - n{n + a) - tH',] + [(2n + a- t)H'^ + H^f. (4.11) 

Proof: First we express r„(t) and /3„(t) in terms of Hn and its derivatives. From (3.24) 
and (4.2) we have 

- Hn = -trn + Pn- n{n + a) 

= -p^{n,t) -n{n + a). (4.12) 

Taking a derivative of (4.2) with respect to t and recalhng (4.3) we have 

r„ = H'^, (4.13) 

and with the first equahty of (4.12) and (4.13), we find, 

pn = tK-Hn + n{n + a). (4.14) 

Now a derivative of (4.14) with respect to t and (Ti) gives 

(tH'J' - H'^ = tH!^ 

= Pn= (-Rji-l — Rn)Pn 
2 

= ^-PnRn- (4.15) 

Here we have made use of (3.13) to arrive at the last equahty. Therefore we have a quadratic 
equation in i?„; 

^-(3^Rn = tH:. (4.16) 

There is another quadratic equation in Rn which is a restatement of (3.16); 

^ + (3nRn = f3n- (2n + a)r„ - n{n + a). (4.17) 

Now we solve for i?„ and 1/-R„ from (4.16) and (4.17) and find 

2r^ 

—^ = pn - {'2n + a)rn - n{n + a) + tH'^ 

2(3nRn = Pn - {2n + a)rn - n{n + a) - tH'^. 
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The equation (4.10) follows from the product of the above two equations, 

Apnrl = [(3n - (2n + a)r„ - n{n + a)]^ - {tHy , 
and (4.13) and (4.14). D 

Incidentally Rn has two alternative representations, 

2[H^ - n{n + a) - tH'^] ^^''''^ 



1 _ tH'^ -{2n + a- t)H'„ - H^ 



(4.19) 



The "discrete" structure inherited from the recurrence relations (2.7), induces a discrete 
analog of the a form, namely, a non-linear second order difference equation in n satisfied 
by Hn for a fixed t ; we believe such a discrete form is new and may have been missed 
in previous similar studies perhaps because the recurrence relations were not sufficiently 
exploited. We note here that our derivation of (4.11) bypasses a third order equation and 
without having to identify a first integral which reduces the order by one. 

We note also that equation (4.11) was first discovered by Tracy and Widom in |27j (which 
in our problem corresponds to A = and B = 1) and just as was done in their paper for 
the Hermite case one can also rescale to obtain the Painleve III equation corresponding to 
the Bessel kernel or "hard edge scaling". We change variables t —>■ s/An, Hn —>■ cr use (4.11) 
and keep only the highest order terms to obtain 

{sa"Y = 4 aia'Y - 4 s(a')' - s{a'Y + aa' + a\a'f. 

Finally, we point out that the above analysis shows that the resolvent kernel used in the 
Tracy- Widom approach can be directly related to the orthogonal polynomials defined on 
(t, oo) . In fact, if we denote R{t,t) as the resolvent kernel defined in [27] then 

tR{t, t) = Hn{t) = -trn - n{n + a) + pn- 
Thus 



hn-lit) K-l{t) 

The term /?„ can also be written using (3.16). In addtion, we have that 



(ti?(t,t)) =r„ 



dt\ ' / K-iit) 

In other words we have found an identity for the resolvent kernel in terms of the values at 
the end points of the normalized orthogonal polynomials. 
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Theorem 7 The auxiliary quantities Rn and r„ are expressed in terms Hn and Hn±i as 
follows: 

t Rn = Hn — Hn+l (4.20) 

_ [Hn - n{n + a)] (t + Hn+i - iJn-i) + tn{n + a) fA 21^ 

t + Hn+i - Hn~\ - 2n - a 

The discrete analog of the a form satisfied by Hn results from the substitution of (4-20) and 
(4.21) into 

{t r„)' = [n{n + a) + t r„ - iJ„][(t i?„)' + t Rn{Hn+i + i/„_i - 2Hn)]. (4.22) 

Proof: Taking a first order difference on the second equality of (4.2) together with (2.10) 
and (3.9) implies (4.20). 
We re-write (3.24) gives 

(3n = n{n + a)+trn- Hn- (4.23) 

We will now find another equation expressing /3„ in terms r„, i?„, if„, Hn±i. Taking a 
first order difference on (4.20) gives, 

t{Rn — Rn-l) = '^Hn — Hn+l — Hn-l- 

Now multiply the above equation by i?„ and make use of (3.13) we find 

tRn JT' ~ V'-^n ~ Hn+l ~ Hn^l)Rni 

Pn 

and therefore 

1 ti?„ — {2Hn — Hn+l — Hn-l)Rn 



Pn tr 



(4.24) 



Therefore the product of (4.23) and (4.24) implies (4.21), which leaves us the job of finding 
a further expression of r„ in terms of Hn and Hn±i- For this purpose we rewrite (3.17) as 

/3„(1 - Rn- Rn-l) = (2n + a)r„ + n{n + a). 

Now substitute /9„ given in (4.23) into the above resulting a linear equation in r„; 

r„[(t - tRn - tRn-i) -2n- a] = [Hn - n{n + a)]{l - Rn - Rn-i) + n(n + a). 

With tRn as in (4.20) we have (4.21). 

We summarize our results in the next theorem 



Theorem 8 Let Dn{t) he the Hankel determinant associated with the Laguerre weight per- 
turbed by a jump factor, and 

Hn{t) ■.= t-\nDnit). 
Then the recurrence coefficients are 

Mt)-(2n + a + l) = ^ (4-25) 

/5„(t)-n(n + a) = tHl^ - H^, (4.26) 

where 2n+l+a and n{n+a) are the "unperturbed" recurrence coefficients and Hn satisfies 
a non-linear differential equation in Jimbo-Miwa-Okamoto a form, 

{tH'^^f = m'n?[Hn - n{n + a) - tH'^] + [{2n + a- t)H', + H^f . 

For the same Hn-, the recurrence coefficients are 

an{t)-{2n + a + l) = H^ - H^+i (4.27) 

Pn{t)-nin + a) = Hni2n + a)- nin + a)iH^,, - H^-,) ^ 

t + Hn+i — Hn-i — 2n — a 

where Hn satisfies the discrete a form of a non-linear difference equation, 

^ 2 
[Hn - n{n + a)]{t + Hn+i - Hn-i) +tn{n + a) 



t + Hn+i — Hn^i — 2n — a 

(2n + a) [Hn — n{n + a)] + tn{n + a) 
t + Hn+i - Hn-1 -2n-a 



■{Hn — Hn+l){Hn-l — Hn). (4.29) 



Note that since Q;„(t) and /3„(t) have two ahernative representations, iJ„(t) satisfies two 
more differential-difference equations, (4.25) = (4.27) and (4.26) = (4.28). 

We end this paper with a discussion on the relationship between our Py and the difference 
equations (3.19) and (3.20). We would like to thank the second referee for supplying us the 
background material part of which is reproduced here. 

The fifth Painleve equation Pyia, b,c,d = —1/2) : 

y = \7r + 7 ]yy) - T + 1^ \ay + -\ +c- + d- 



2y y-lj^''' t t2 \" yj t 1/ - 1 ' dt 

is equivalent to the Hamiltonian system Tiy '■ 

op oq 
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and 



with the time-dependent Hamihonian H = \-\{p,q,t) : 

tH = p{p + t)q{q - 1) + a2qt - a^pq - aip{q - 1), 

where 

a = — , b = — — , c = ao-a2, a = --, ao := 1 - ai - 02 - "3 

1 

The Hamiltonian structure was studied in p^ and the r— function is defined such that 

4lnr = H. 

dt 

The extended affine Weyl group ^^(Ag ') =< sq, Si, S2, S3, vr > of the Weyl group type Ag 
acts as bi-rational symmetries on Py and induces Backlund transformations on the solutions 
of Py- Here the s[ s and vr are the generators. See [22] for the study of Weyl group actions 
on Pv . 

For example, the action of sq : 

so{ao, ai, a2, 03} = {-ao, ai + ao, "2, ^3 + "o}, 

leaves Py or the Hamiltonian system Tiy invariant. We refer the readers to [221 [2H] for 
information on Weyl group actions and [(4.3),|[28j] which lists the bi-rational transformations. 
To proceed further, consider a parallel transformation / = (s2S37r)^ G H^(^3 ) : 

I : a = {ao, ai, a2, 03) 1 — ^ a + (1, 0, —1, 0). 

From a direct computation, we may verify that, the variables q and r '■= pq{q — I) satisfies 
the following system of difference equations: 

/(r) + r = q{a2 — a^ + t — tq) — ai (4.30) 



1 _ \ / 1 _ ;L^ ^ {r -a2){r-a2-a3) 
Q J \l~^{(l) J r{r + ai) 



(4.31) 



these seems to the d — Pni of [205,|i23|] in disguise. 
In terms of H our auxiliary parameter r reads 



±m. (4.32) 
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In our problem we have Py(0, —a^/2, 2?7, + 1 + a, —1/2), which imphes 

(ao, ai, 0L2, a-i) = (1 + n, 0, —n — a, a). 
If 

/"(r) =: r„, 
for n G {0, 1,2, ...}, then a direct computation shows that (4.33) and (4.34) 



are 



r„+i + r„ = Rn{-a-2n~l+t-tRn), (4.33) 



1 .\ / 1 ,\ _ {rn + n + a){rn + n) 



i?n y V-Rn-1 J rl 



(4.34) 



are equivalent to (3.19) and (3.20) respectively. 

We should like to mention here that (3.19) and (3.20) and other equations are derived 
entirely from orthogonality and the immediate consequence — the recurrence relations. 

In view of (3.32) we see that the logarithmic derivative of the generating function 
G(n, t) = Dn{t) is the r— function of our Py- We end this paper with the final remark: the 
equation (4.4) is essentially the same as the Toda equation among a r— sequence discovered 
by Okamoto |22j . 
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